3.643 \(\int \frac {1}{(d f+e f x)^2 (a+b (d+e x)^2+c (d+e x)^4)} \, dx\)

Optimal. Leaf size=204 \[ -\frac {\sqrt {c} \left (\frac {b}{\sqrt {b^2-4 a c}}+1\right ) \tan ^{-1}\left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {b-\sqrt {b^2-4 a c}}}\right )}{\sqrt {2} a e f^2 \sqrt {b-\sqrt {b^2-4 a c}}}-\frac {\sqrt {c} \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right ) \tan ^{-1}\left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {\sqrt {b^2-4 a c}+b}}\right )}{\sqrt {2} a e f^2 \sqrt {\sqrt {b^2-4 a c}+b}}-\frac {1}{a e f^2 (d+e x)} \]

[Out]

-1/a/e/f^2/(e*x+d)-1/2*arctan((e*x+d)*2^(1/2)*c^(1/2)/(b-(-4*a*c+b^2)^(1/2))^(1/2))*c^(1/2)*(1+b/(-4*a*c+b^2)^
(1/2))/a/e/f^2*2^(1/2)/(b-(-4*a*c+b^2)^(1/2))^(1/2)-1/2*arctan((e*x+d)*2^(1/2)*c^(1/2)/(b+(-4*a*c+b^2)^(1/2))^
(1/2))*c^(1/2)*(1-b/(-4*a*c+b^2)^(1/2))/a/e/f^2*2^(1/2)/(b+(-4*a*c+b^2)^(1/2))^(1/2)

________________________________________________________________________________________

Rubi [A]  time = 0.27, antiderivative size = 204, normalized size of antiderivative = 1.00, number of steps used = 5, number of rules used = 4, integrand size = 33, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.121, Rules used = {1142, 1123, 1166, 205} \[ -\frac {\sqrt {c} \left (\frac {b}{\sqrt {b^2-4 a c}}+1\right ) \tan ^{-1}\left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {b-\sqrt {b^2-4 a c}}}\right )}{\sqrt {2} a e f^2 \sqrt {b-\sqrt {b^2-4 a c}}}-\frac {\sqrt {c} \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right ) \tan ^{-1}\left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {\sqrt {b^2-4 a c}+b}}\right )}{\sqrt {2} a e f^2 \sqrt {\sqrt {b^2-4 a c}+b}}-\frac {1}{a e f^2 (d+e x)} \]

Antiderivative was successfully verified.

[In]

Int[1/((d*f + e*f*x)^2*(a + b*(d + e*x)^2 + c*(d + e*x)^4)),x]

[Out]

-(1/(a*e*f^2*(d + e*x))) - (Sqrt[c]*(1 + b/Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*(d + e*x))/Sqrt[b - Sqrt
[b^2 - 4*a*c]]])/(Sqrt[2]*a*Sqrt[b - Sqrt[b^2 - 4*a*c]]*e*f^2) - (Sqrt[c]*(1 - b/Sqrt[b^2 - 4*a*c])*ArcTan[(Sq
rt[2]*Sqrt[c]*(d + e*x))/Sqrt[b + Sqrt[b^2 - 4*a*c]]])/(Sqrt[2]*a*Sqrt[b + Sqrt[b^2 - 4*a*c]]*e*f^2)

Rule 205

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]*ArcTan[x/Rt[a/b, 2]])/a, x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]

Rule 1123

Int[((d_.)*(x_))^(m_)*((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4)^(p_), x_Symbol] :> Simp[((d*x)^(m + 1)*(a + b*x^2 +
 c*x^4)^(p + 1))/(a*d*(m + 1)), x] - Dist[1/(a*d^2*(m + 1)), Int[(d*x)^(m + 2)*(b*(m + 2*p + 3) + c*(m + 4*p +
 5)*x^2)*(a + b*x^2 + c*x^4)^p, x], x] /; FreeQ[{a, b, c, d, p}, x] && NeQ[b^2 - 4*a*c, 0] && LtQ[m, -1] && In
tegerQ[2*p] && (IntegerQ[p] || IntegerQ[m])

Rule 1142

Int[(u_)^(m_.)*((a_.) + (b_.)*(v_)^2 + (c_.)*(v_)^4)^(p_.), x_Symbol] :> Dist[u^m/(Coefficient[v, x, 1]*v^m),
Subst[Int[x^m*(a + b*x^2 + c*x^(2*2))^p, x], x, v], x] /; FreeQ[{a, b, c, m, p}, x] && LinearPairQ[u, v, x]

Rule 1166

Int[((d_) + (e_.)*(x_)^2)/((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[b^2 - 4*a*c, 2]}, Di
st[e/2 + (2*c*d - b*e)/(2*q), Int[1/(b/2 - q/2 + c*x^2), x], x] + Dist[e/2 - (2*c*d - b*e)/(2*q), Int[1/(b/2 +
 q/2 + c*x^2), x], x]] /; FreeQ[{a, b, c, d, e}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - a*e^2, 0] && PosQ[b^
2 - 4*a*c]

Rubi steps

\begin {align*} \int \frac {1}{(d f+e f x)^2 \left (a+b (d+e x)^2+c (d+e x)^4\right )} \, dx &=\frac {\operatorname {Subst}\left (\int \frac {1}{x^2 \left (a+b x^2+c x^4\right )} \, dx,x,d+e x\right )}{e f^2}\\ &=-\frac {1}{a e f^2 (d+e x)}+\frac {\operatorname {Subst}\left (\int \frac {-b-c x^2}{a+b x^2+c x^4} \, dx,x,d+e x\right )}{a e f^2}\\ &=-\frac {1}{a e f^2 (d+e x)}-\frac {\left (c \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right )\right ) \operatorname {Subst}\left (\int \frac {1}{\frac {b}{2}+\frac {1}{2} \sqrt {b^2-4 a c}+c x^2} \, dx,x,d+e x\right )}{2 a e f^2}-\frac {\left (c \left (1+\frac {b}{\sqrt {b^2-4 a c}}\right )\right ) \operatorname {Subst}\left (\int \frac {1}{\frac {b}{2}-\frac {1}{2} \sqrt {b^2-4 a c}+c x^2} \, dx,x,d+e x\right )}{2 a e f^2}\\ &=-\frac {1}{a e f^2 (d+e x)}-\frac {\sqrt {c} \left (1+\frac {b}{\sqrt {b^2-4 a c}}\right ) \tan ^{-1}\left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {b-\sqrt {b^2-4 a c}}}\right )}{\sqrt {2} a \sqrt {b-\sqrt {b^2-4 a c}} e f^2}-\frac {\sqrt {c} \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right ) \tan ^{-1}\left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {b+\sqrt {b^2-4 a c}}}\right )}{\sqrt {2} a \sqrt {b+\sqrt {b^2-4 a c}} e f^2}\\ \end {align*}

________________________________________________________________________________________

Mathematica [A]  time = 0.34, size = 209, normalized size = 1.02 \[ -\frac {\frac {\sqrt {2} \sqrt {c} \left (\sqrt {b^2-4 a c}+b\right ) \tan ^{-1}\left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {b-\sqrt {b^2-4 a c}}}\right )}{\sqrt {b^2-4 a c} \sqrt {b-\sqrt {b^2-4 a c}}}+\frac {\sqrt {2} \sqrt {c} \left (\sqrt {b^2-4 a c}-b\right ) \tan ^{-1}\left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {\sqrt {b^2-4 a c}+b}}\right )}{\sqrt {b^2-4 a c} \sqrt {\sqrt {b^2-4 a c}+b}}+\frac {2}{d+e x}}{2 a e f^2} \]

Antiderivative was successfully verified.

[In]

Integrate[1/((d*f + e*f*x)^2*(a + b*(d + e*x)^2 + c*(d + e*x)^4)),x]

[Out]

-1/2*(2/(d + e*x) + (Sqrt[2]*Sqrt[c]*(b + Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*(d + e*x))/Sqrt[b - Sqrt[
b^2 - 4*a*c]]])/(Sqrt[b^2 - 4*a*c]*Sqrt[b - Sqrt[b^2 - 4*a*c]]) + (Sqrt[2]*Sqrt[c]*(-b + Sqrt[b^2 - 4*a*c])*Ar
cTan[(Sqrt[2]*Sqrt[c]*(d + e*x))/Sqrt[b + Sqrt[b^2 - 4*a*c]]])/(Sqrt[b^2 - 4*a*c]*Sqrt[b + Sqrt[b^2 - 4*a*c]])
)/(a*e*f^2)

________________________________________________________________________________________

fricas [B]  time = 0.94, size = 1477, normalized size = 7.24 \[ \text {result too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(e*f*x+d*f)^2/(a+b*(e*x+d)^2+c*(e*x+d)^4),x, algorithm="fricas")

[Out]

1/2*(sqrt(1/2)*(a*e^2*f^2*x + a*d*e*f^2)*sqrt(-((a^3*b^2 - 4*a^4*c)*e^2*f^4*sqrt((b^4 - 2*a*b^2*c + a^2*c^2)/(
(a^6*b^2 - 4*a^7*c)*e^4*f^8)) + b^3 - 3*a*b*c)/((a^3*b^2 - 4*a^4*c)*e^2*f^4))*log(-2*(b^2*c^2 - a*c^3)*e*x - 2
*(b^2*c^2 - a*c^3)*d + sqrt(1/2)*((a^3*b^4 - 6*a^4*b^2*c + 8*a^5*c^2)*e^3*f^6*sqrt((b^4 - 2*a*b^2*c + a^2*c^2)
/((a^6*b^2 - 4*a^7*c)*e^4*f^8)) - (b^5 - 5*a*b^3*c + 4*a^2*b*c^2)*e*f^2)*sqrt(-((a^3*b^2 - 4*a^4*c)*e^2*f^4*sq
rt((b^4 - 2*a*b^2*c + a^2*c^2)/((a^6*b^2 - 4*a^7*c)*e^4*f^8)) + b^3 - 3*a*b*c)/((a^3*b^2 - 4*a^4*c)*e^2*f^4)))
 - sqrt(1/2)*(a*e^2*f^2*x + a*d*e*f^2)*sqrt(-((a^3*b^2 - 4*a^4*c)*e^2*f^4*sqrt((b^4 - 2*a*b^2*c + a^2*c^2)/((a
^6*b^2 - 4*a^7*c)*e^4*f^8)) + b^3 - 3*a*b*c)/((a^3*b^2 - 4*a^4*c)*e^2*f^4))*log(-2*(b^2*c^2 - a*c^3)*e*x - 2*(
b^2*c^2 - a*c^3)*d - sqrt(1/2)*((a^3*b^4 - 6*a^4*b^2*c + 8*a^5*c^2)*e^3*f^6*sqrt((b^4 - 2*a*b^2*c + a^2*c^2)/(
(a^6*b^2 - 4*a^7*c)*e^4*f^8)) - (b^5 - 5*a*b^3*c + 4*a^2*b*c^2)*e*f^2)*sqrt(-((a^3*b^2 - 4*a^4*c)*e^2*f^4*sqrt
((b^4 - 2*a*b^2*c + a^2*c^2)/((a^6*b^2 - 4*a^7*c)*e^4*f^8)) + b^3 - 3*a*b*c)/((a^3*b^2 - 4*a^4*c)*e^2*f^4))) -
 sqrt(1/2)*(a*e^2*f^2*x + a*d*e*f^2)*sqrt(((a^3*b^2 - 4*a^4*c)*e^2*f^4*sqrt((b^4 - 2*a*b^2*c + a^2*c^2)/((a^6*
b^2 - 4*a^7*c)*e^4*f^8)) - b^3 + 3*a*b*c)/((a^3*b^2 - 4*a^4*c)*e^2*f^4))*log(-2*(b^2*c^2 - a*c^3)*e*x - 2*(b^2
*c^2 - a*c^3)*d + sqrt(1/2)*((a^3*b^4 - 6*a^4*b^2*c + 8*a^5*c^2)*e^3*f^6*sqrt((b^4 - 2*a*b^2*c + a^2*c^2)/((a^
6*b^2 - 4*a^7*c)*e^4*f^8)) + (b^5 - 5*a*b^3*c + 4*a^2*b*c^2)*e*f^2)*sqrt(((a^3*b^2 - 4*a^4*c)*e^2*f^4*sqrt((b^
4 - 2*a*b^2*c + a^2*c^2)/((a^6*b^2 - 4*a^7*c)*e^4*f^8)) - b^3 + 3*a*b*c)/((a^3*b^2 - 4*a^4*c)*e^2*f^4))) + sqr
t(1/2)*(a*e^2*f^2*x + a*d*e*f^2)*sqrt(((a^3*b^2 - 4*a^4*c)*e^2*f^4*sqrt((b^4 - 2*a*b^2*c + a^2*c^2)/((a^6*b^2
- 4*a^7*c)*e^4*f^8)) - b^3 + 3*a*b*c)/((a^3*b^2 - 4*a^4*c)*e^2*f^4))*log(-2*(b^2*c^2 - a*c^3)*e*x - 2*(b^2*c^2
 - a*c^3)*d - sqrt(1/2)*((a^3*b^4 - 6*a^4*b^2*c + 8*a^5*c^2)*e^3*f^6*sqrt((b^4 - 2*a*b^2*c + a^2*c^2)/((a^6*b^
2 - 4*a^7*c)*e^4*f^8)) + (b^5 - 5*a*b^3*c + 4*a^2*b*c^2)*e*f^2)*sqrt(((a^3*b^2 - 4*a^4*c)*e^2*f^4*sqrt((b^4 -
2*a*b^2*c + a^2*c^2)/((a^6*b^2 - 4*a^7*c)*e^4*f^8)) - b^3 + 3*a*b*c)/((a^3*b^2 - 4*a^4*c)*e^2*f^4))) - 2)/(a*e
^2*f^2*x + a*d*e*f^2)

________________________________________________________________________________________

giac [F(-2)]  time = 0.00, size = 0, normalized size = 0.00 \[ \text {Exception raised: TypeError} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(e*f*x+d*f)^2/(a+b*(e*x+d)^2+c*(e*x+d)^4),x, algorithm="giac")

[Out]

Exception raised: TypeError >> An error occurred running a Giac command:INPUT:sage2:=int(sage0,x):;OUTPUT:Erro
r index.cc index_gcd Error: Bad Argument ValueError index.cc index_gcd Error: Bad Argument ValueDone

________________________________________________________________________________________

maple [C]  time = 0.01, size = 174, normalized size = 0.85 \[ \frac {\left (-\RootOf \left (\textit {\_Z}^{4} c \,e^{4}+4 \textit {\_Z}^{3} c d \,e^{3}+c \,d^{4}+b \,d^{2}+\left (6 c \,d^{2} e^{2}+b \,e^{2}\right ) \textit {\_Z}^{2}+\left (4 c \,d^{3} e +2 d e b \right ) \textit {\_Z} +a \right )^{2} c \,e^{2}-2 \RootOf \left (\textit {\_Z}^{4} c \,e^{4}+4 \textit {\_Z}^{3} c d \,e^{3}+c \,d^{4}+b \,d^{2}+\left (6 c \,d^{2} e^{2}+b \,e^{2}\right ) \textit {\_Z}^{2}+\left (4 c \,d^{3} e +2 d e b \right ) \textit {\_Z} +a \right ) c d e -c \,d^{2}-b \right ) \ln \left (-\RootOf \left (\textit {\_Z}^{4} c \,e^{4}+4 \textit {\_Z}^{3} c d \,e^{3}+c \,d^{4}+b \,d^{2}+\left (6 c \,d^{2} e^{2}+b \,e^{2}\right ) \textit {\_Z}^{2}+\left (4 c \,d^{3} e +2 d e b \right ) \textit {\_Z} +a \right )+x \right )}{2 a e \,f^{2} \left (2 c \,e^{3} \RootOf \left (\textit {\_Z}^{4} c \,e^{4}+4 \textit {\_Z}^{3} c d \,e^{3}+c \,d^{4}+b \,d^{2}+\left (6 c \,d^{2} e^{2}+b \,e^{2}\right ) \textit {\_Z}^{2}+\left (4 c \,d^{3} e +2 d e b \right ) \textit {\_Z} +a \right )^{3}+6 c d \,e^{2} \RootOf \left (\textit {\_Z}^{4} c \,e^{4}+4 \textit {\_Z}^{3} c d \,e^{3}+c \,d^{4}+b \,d^{2}+\left (6 c \,d^{2} e^{2}+b \,e^{2}\right ) \textit {\_Z}^{2}+\left (4 c \,d^{3} e +2 d e b \right ) \textit {\_Z} +a \right )^{2}+6 e c \,d^{2} \RootOf \left (\textit {\_Z}^{4} c \,e^{4}+4 \textit {\_Z}^{3} c d \,e^{3}+c \,d^{4}+b \,d^{2}+\left (6 c \,d^{2} e^{2}+b \,e^{2}\right ) \textit {\_Z}^{2}+\left (4 c \,d^{3} e +2 d e b \right ) \textit {\_Z} +a \right )+2 c \,d^{3}+b e \RootOf \left (\textit {\_Z}^{4} c \,e^{4}+4 \textit {\_Z}^{3} c d \,e^{3}+c \,d^{4}+b \,d^{2}+\left (6 c \,d^{2} e^{2}+b \,e^{2}\right ) \textit {\_Z}^{2}+\left (4 c \,d^{3} e +2 d e b \right ) \textit {\_Z} +a \right )+b d \right )}-\frac {1}{\left (e x +d \right ) a e \,f^{2}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(1/(e*f*x+d*f)^2/(a+b*(e*x+d)^2+c*(e*x+d)^4),x)

[Out]

1/2/f^2/a/e*sum((-_R^2*c*e^2-2*_R*c*d*e-c*d^2-b)/(2*_R^3*c*e^3+6*_R^2*c*d*e^2+6*_R*c*d^2*e+2*c*d^3+_R*b*e+b*d)
*ln(-_R+x),_R=RootOf(_Z^4*c*e^4+4*_Z^3*c*d*e^3+c*d^4+b*d^2+(6*c*d^2*e^2+b*e^2)*_Z^2+(4*c*d^3*e+2*b*d*e)*_Z+a))
-1/a/e/f^2/(e*x+d)

________________________________________________________________________________________

maxima [F(-1)]  time = 0.00, size = 0, normalized size = 0.00 \[ \text {Timed out} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(e*f*x+d*f)^2/(a+b*(e*x+d)^2+c*(e*x+d)^4),x, algorithm="maxima")

[Out]

Timed out

________________________________________________________________________________________

mupad [B]  time = 3.87, size = 4339, normalized size = 21.27 \[ \text {result too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(1/((d*f + e*f*x)^2*(a + b*(d + e*x)^2 + c*(d + e*x)^4)),x)

[Out]

- atan(((-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^(1/2))/(8*(a
^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2)*(x*(4*a^4*c^4*e^12*f^6 - 2*a^3*b^2*c^3*e^12
*f^6) - ((x*(8*a^5*b^3*c^2*e^14*f^10 - 32*a^6*b*c^3*e^14*f^10) - 32*a^6*b*c^3*d*e^13*f^10 + 8*a^5*b^3*c^2*d*e^
13*f^10)*(-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^(1/2))/(8*(
a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) - 4*a^4*b^3*c^2*e^12*f^8 + 16*a^5*b*c^3*e^
12*f^8)*(-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^(1/2))/(8*(a
^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*c^4*d*e^11*f^6 - 2*a^3*b^2*c^3*d*e^
11*f^6)*1i + (-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^(1/2))/
(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2)*(x*(4*a^4*c^4*e^12*f^6 - 2*a^3*b^2*c^3
*e^12*f^6) - ((x*(8*a^5*b^3*c^2*e^14*f^10 - 32*a^6*b*c^3*e^14*f^10) - 32*a^6*b*c^3*d*e^13*f^10 + 8*a^5*b^3*c^2
*d*e^13*f^10)*(-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^(1/2))
/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*b^3*c^2*e^12*f^8 - 16*a^5*b*c
^3*e^12*f^8)*(-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^(1/2))/
(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*c^4*d*e^11*f^6 - 2*a^3*b^2*c^3
*d*e^11*f^6)*1i)/((-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^(1
/2))/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2)*(x*(4*a^4*c^4*e^12*f^6 - 2*a^3*b^
2*c^3*e^12*f^6) - ((x*(8*a^5*b^3*c^2*e^14*f^10 - 32*a^6*b*c^3*e^14*f^10) - 32*a^6*b*c^3*d*e^13*f^10 + 8*a^5*b^
3*c^2*d*e^13*f^10)*(-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^(
1/2))/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*b^3*c^2*e^12*f^8 - 16*a^
5*b*c^3*e^12*f^8)*(-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^(1
/2))/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*c^4*d*e^11*f^6 - 2*a^3*b^
2*c^3*d*e^11*f^6) - (-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^
(1/2))/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2)*(x*(4*a^4*c^4*e^12*f^6 - 2*a^3*
b^2*c^3*e^12*f^6) - ((x*(8*a^5*b^3*c^2*e^14*f^10 - 32*a^6*b*c^3*e^14*f^10) - 32*a^6*b*c^3*d*e^13*f^10 + 8*a^5*
b^3*c^2*d*e^13*f^10)*(-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)
^(1/2))/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) - 4*a^4*b^3*c^2*e^12*f^8 + 16*
a^5*b*c^3*e^12*f^8)*(-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^
(1/2))/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*c^4*d*e^11*f^6 - 2*a^3*
b^2*c^3*d*e^11*f^6) + 2*a^3*c^4*e^10*f^4))*(-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c -
a*c*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2)*2i - ata
n(((-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^3*b^
4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2)*(x*(4*a^4*c^4*e^12*f^6 - 2*a^3*b^2*c^3*e^12*f^6)
 - ((x*(8*a^5*b^3*c^2*e^14*f^10 - 32*a^6*b*c^3*e^14*f^10) - 32*a^6*b*c^3*d*e^13*f^10 + 8*a^5*b^3*c^2*d*e^13*f^
10)*(-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^3*b
^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) - 4*a^4*b^3*c^2*e^12*f^8 + 16*a^5*b*c^3*e^12*f^
8)*(-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^3*b^
4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*c^4*d*e^11*f^6 - 2*a^3*b^2*c^3*d*e^11*f^
6)*1i + (-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2))/(8*(a
^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2)*(x*(4*a^4*c^4*e^12*f^6 - 2*a^3*b^2*c^3*e^12
*f^6) - ((x*(8*a^5*b^3*c^2*e^14*f^10 - 32*a^6*b*c^3*e^14*f^10) - 32*a^6*b*c^3*d*e^13*f^10 + 8*a^5*b^3*c^2*d*e^
13*f^10)*(-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2))/(8*(
a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*b^3*c^2*e^12*f^8 - 16*a^5*b*c^3*e^
12*f^8)*(-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2))/(8*(a
^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*c^4*d*e^11*f^6 - 2*a^3*b^2*c^3*d*e^
11*f^6)*1i)/((-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2))/
(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2)*(x*(4*a^4*c^4*e^12*f^6 - 2*a^3*b^2*c^3
*e^12*f^6) - ((x*(8*a^5*b^3*c^2*e^14*f^10 - 32*a^6*b*c^3*e^14*f^10) - 32*a^6*b*c^3*d*e^13*f^10 + 8*a^5*b^3*c^2
*d*e^13*f^10)*(-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2))
/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*b^3*c^2*e^12*f^8 - 16*a^5*b*c
^3*e^12*f^8)*(-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2))/
(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*c^4*d*e^11*f^6 - 2*a^3*b^2*c^3
*d*e^11*f^6) - (-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2)
)/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2)*(x*(4*a^4*c^4*e^12*f^6 - 2*a^3*b^2*c
^3*e^12*f^6) - ((x*(8*a^5*b^3*c^2*e^14*f^10 - 32*a^6*b*c^3*e^14*f^10) - 32*a^6*b*c^3*d*e^13*f^10 + 8*a^5*b^3*c
^2*d*e^13*f^10)*(-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2
))/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) - 4*a^4*b^3*c^2*e^12*f^8 + 16*a^5*b
*c^3*e^12*f^8)*(-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2)
)/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*c^4*d*e^11*f^6 - 2*a^3*b^2*c
^3*d*e^11*f^6) + 2*a^3*c^4*e^10*f^4))*(-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(
-(4*a*c - b^2)^3)^(1/2))/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2)*2i - 1/(a*e*(
d*f^2 + e*f^2*x))

________________________________________________________________________________________

sympy [A]  time = 4.90, size = 258, normalized size = 1.26 \[ \operatorname {RootSum} {\left (t^{4} \left (256 a^{5} c^{2} e^{4} f^{8} - 128 a^{4} b^{2} c e^{4} f^{8} + 16 a^{3} b^{4} e^{4} f^{8}\right ) + t^{2} \left (48 a^{2} b c^{2} e^{2} f^{4} - 28 a b^{3} c e^{2} f^{4} + 4 b^{5} e^{2} f^{4}\right ) + c^{3}, \left (t \mapsto t \log {\left (x + \frac {- 64 t^{3} a^{5} c^{2} e^{3} f^{6} + 48 t^{3} a^{4} b^{2} c e^{3} f^{6} - 8 t^{3} a^{3} b^{4} e^{3} f^{6} - 10 t a^{2} b c^{2} e f^{2} + 10 t a b^{3} c e f^{2} - 2 t b^{5} e f^{2} + a c^{3} d - b^{2} c^{2} d}{a c^{3} e - b^{2} c^{2} e} \right )} \right )\right )} - \frac {1}{a d e f^{2} + a e^{2} f^{2} x} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(1/(e*f*x+d*f)**2/(a+b*(e*x+d)**2+c*(e*x+d)**4),x)

[Out]

RootSum(_t**4*(256*a**5*c**2*e**4*f**8 - 128*a**4*b**2*c*e**4*f**8 + 16*a**3*b**4*e**4*f**8) + _t**2*(48*a**2*
b*c**2*e**2*f**4 - 28*a*b**3*c*e**2*f**4 + 4*b**5*e**2*f**4) + c**3, Lambda(_t, _t*log(x + (-64*_t**3*a**5*c**
2*e**3*f**6 + 48*_t**3*a**4*b**2*c*e**3*f**6 - 8*_t**3*a**3*b**4*e**3*f**6 - 10*_t*a**2*b*c**2*e*f**2 + 10*_t*
a*b**3*c*e*f**2 - 2*_t*b**5*e*f**2 + a*c**3*d - b**2*c**2*d)/(a*c**3*e - b**2*c**2*e)))) - 1/(a*d*e*f**2 + a*e
**2*f**2*x)

________________________________________________________________________________________